Island Distance in One— Dimensional Epitaxial Growth 
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The typical island distance I in submonlayer epitaxial growth depends on the growth conditions via 
an exponent 7. This exponent is known to depend on the substrate dimensionality, the dimension 
of the islands, and the size i* of the critical nucleus for island formation. In this paper we study 
the dependence of 7 on i* in one-dimensional epitaxial growth. We derive that 7 = z*/(2i* + 3) for 
i* > 2 and confirm this result by computer simulations. 

PACS numbers: 81.15-z, 81.10.Bk, 05.50.+q 

I. INTRODUCTION 

Molecular beam epitaxy (MBE) is of outstanding im- 
portance for the fabrication of microelectronic devices. 
The highly controllable experimental conditions make 
this technique an ideal tool for crystal growth. The thick- 
ness of the deposited layers can be controlled down to the 
atomic scale while growth proceeds layerwise. 

Typically the growth of multiple layers is of primary 
interest in most experimental and theoretical studies of 
MBE. However, deeper insight into the physics of multi- 
layer growth requires a comprehensive understanding of 
the submonolayer regime (see Ref. for instance). The 
most important length scale in submonolayer growth is 
the island distance or diffusion length I. The classical re- 
sult B for this length is that it depends on the diffusion 
constant D of adatoms and the deposition rate F like 



(3) 
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The analysis of rate equations yields that the exponent 
7 depends on the substrate dimension d, on the size i* of 
the critical nucleus and the fractal dimension df of the 
islands MM: 
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d + df 



(2) 



The critical nucleus is defined such that all islands con- 
sisting of i* or less atoms can shrink by emitting adatoms, 
while an island of size i* + 1 is stable. Eq. (Q) applies 
if island diffusion and adatom desorption are not impor- 
tant. 

Eq. (||) has been demonstrated to hold in d = 2 for 
i* = 1,2,3,4, see Refs. In this paper we study the 

dependence (Q) in surface dimension d = 1 for different 
values of the critical nucleus i*. In this case, df — 1. So 
far, only the case i* = 1 ^|J7) has been studied in one 
dimension, to our knowledge, and it is in agreement with 
(0). Here we show that (0) is no longer valid for i* > 2, 
where 



Our motivation for studying the one-dimensional case 
is the following. Much theoretical work is dedicated to 
the study of basic mechanisms of epitaxial growth. For 
computational reasons, the numerical work in this area 
is often restricted to models in surface dimension d = 1 . 
As mentioned above, such studies rely on a good under- 
standing of the physics in the submonolayer regime. 

Also, quasi one-dimensional MBE may be realized 
experimentally on surfaces with extremely anisotropic 
diffusion constant, like the 2xl-dimer-reconstructed 
Si(OOl) surface (see, e.g., Ref. ||). There, the esti- 
mated ratio of the diffusion constants along rows and 
perpendicular to them is of the order of 1000:1 in typi- 
cal experiments 0] . A computer simulation study of this 
growth process leads to the estimate that for a deposi- 
tion rate of F = 1/600 ML/s and at temperatures below 
T = 450 K, this anisotropy should be pronounced enough 
to lead to quasi one- dimensional MBE 0. However, 
there is a controversy concerning the relative importance 
of anisotropic diffusion and anisotropic bonding pLpj]]. 
Another experimental example is the reconstructed sur- 
face of Au(100) 0. 

This paper is organized as follows. In section II the 
expression (|^) will be derived. The subsequent sections 
are devoted to the numerical confirmation of this re- 
sult. In section III we employ the coarse-grained MBE 
(CGMBE) model |], which allows for a simple algorith- 
mic implementation of a variable critical nucleus size i*. 
To ensure that the results are independent of the chosen 
algorithm, we reproduce them with a model that is not 
coarse grained (MBE model) in section IV. It is shown 
that the results of the MBE and the CGMBE models 
compare favourably with each other and with the theo- 
retical prediction. The final section V contains a general 
discussion of the results. 
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II. ANALYTICAL CONSIDERATIONS 

Most treatments of the monolayer growth processes 
are based on a mean-field rate equation approach. This 
is appropriate in the most important two-dimensional 
case since the basic kinetic mechanism - the two-particle 
diffusion-controlled aggregation - has an upper critical 
dimension d c = 2 In one dimension, however, clas- 
sic rate equations do not apply and we shall employ a 
modified rate equation approach [Q. 

We denote by A and / densities of adatoms and stable 
immobile islands, respectively. We start with adatoms 
and write 



dA 
dt 
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A 
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Here r is a lifetime, i.e., the time between the deposi- 
tion of an adatom and its incorporation into a stable is- 
land. During the time interval t, the adatom visits \J Dt 
different sites, so the adatom lifetime is determined by 
1\TDt ps 1, and Eq. (0) becomes 



^=F- DI 2 A. 
dt 



(5) 



This differs from the two-dimensional situation where the 
number of different sites visited by an adatom grows lin- 
early with time and the adatom density evolves according 
to dA/dt = F - DIA. Actually, Eq. (|) contains an ad- 
ditional loss term due to the nucleation of new islands. 
However, this loss term is asymptotically subdominant, 
since most adatoms are captured by the existing islands. 

The rate equation for the island density is derived by 
noting that islands are formed upon colliding of i* + 1 
adatoms. Hence 



d -L=DA^\ 
dt 



(6) 



The many-particle diffusion-controlled reaction process 
(i* + l)A — » / has an upper critical dimension d c = 2/i* 
(see e.g. Jl5| and references therein). The fact that the 
upper critical dimension gets smaller when i* increases 
is simple to appreciate. Indeed, many-particle collisions 
occur rarely, so the system is well mixed and the mean- 
field approach should become correct earlier. For i* > 2 
we have d c < 1, so Eq. (^|) applies; the case i* = 2 is 
marginal, so we anticipate logarithmic corrections to the 
mean-field results. For i* = 1 correlations dominate the 
diffusion process in one dimension, so that (^|) has to be 
modified (H) and leads to 7 = 1/4 @. 

Eq.(||) is obtained from (||) and (||) in the follow- 
ing way: From (jsj) it follows that asymptotically A ~ 
(F/D)I~ 2 . Inserting this into (^) and intergrating we 
find 



i*/(2i*+3) 



(Ft) 



l/(2i*+3) 



(7) 



Setting t oc F 1 gives the result (||). 

For i* = 2 there is a logarithmic correction to 

dl _ DA 3 
dt ~ \lnlV 



(8) 



whereas (|^) remains unchanged. Repeating the same 
analysis gives 



2/7 r 



ln 



D 



1/7 
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III. NUMERICAL RESULTS 

The idea of the coarse grained MBE (CGMBE) model, 
introduced in Ref. ||, successfully applied in Refs. H|), 
and studied in detail in Ref. p6| , is to avoid the time- 
consuming simulation of the surface diffusion on the 
atomic scale by resolving it only on a much coarser scale 
Ax. As long as Ax is smaller than the island distance £, 
one still gets correct information about the surface mor- 
phology. 

Specificially the model is implemented as follows: The 
substrate is divided into cells of linear extension Ax, so 
that a monolayer corresponds to Ax atoms per cell. The 
partial filling of a cell at site x is given by a counter 
m(x) = 0, 1, ... , Ax — 1, which means the number of 
atoms on top of completed monolayers in cell x. The 
height h(x) denotes the number of complete monolayers. 
For the complete description of the state of a cell, an edge 
indicator is(x) is introduced. If is(x) = 1, atoms cannot 
leave cell x because they are irreversibly incorporated 
into the island edge present in the cell. This allows for 
a simple implementation of the concept of the critical 
nucleus i*: if upon diffusion or deposition of an atom 
into cell x the counter reaches the value m(x) = i* , is(x) 
is set to one. Consequently, m(x) can only increase from 
then on, until m(x) = Ax, in which case h(x) is increased 
by one and m(x) and is(x) are set to zero. As long as 
is(x) = 0, m(x) can also decrease, when atoms diffuse 
into neighbouring cells. 

The dynamics has to take into account that a single 
adatom move from one cell to a neighbouring one repre- 
sents the diffusion over a distance of Ax lattice sites. This 
happens with frequency v-q = D/(Ax) 2 for all adatoms 
in the system. The deposition of one atom happens with 
frequency v-p = FLAx. Hence the probability for de- 
position during one time step At = (v-p + isp,) 1 of the 
computer simulation is 



(10) 



1 — p is the probability with which each of the adatoms is 
allowed to diffuse into a neighboring cell. If no adatoms 
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are present, = 0. This algorithm speeds the simula- 
tion up by a factor of about (Ace) 2 . 

In the simulations presented in the following, the 
coarse graining length had the constant value Ax = 20. 
The system size was chosen to be LAx = 2 • 10 5 and 
averages were taken over 100 independent runs. 




D/F 

FIG. 1. The island distance as a function of D/F 
for i* = 1,2,3,4. The solid lines represent fits to 
the last three decades of data and have the slopes 
7 = 0.248, 0.293, 0.329, 0.355, respectively. 

In the CGMBE model, nucleation events can be iden- 
tified easily. Thus the island distance can be determined 
by dividing the system size by the total number of nu- 
cleation events in the first layer. Measurements of this 
quantity as a function of D/F for different i* are shown 
in Fig. [|. Note that the asymptotic algebraic behaviour 
£ ~ (D/F) 1 is reached with i as function of D/F having 
a weak negative curvature. 



IV. RESULTS WITHOUT COARSE GRAINING 

In the (non-coarse grained) MBE model atoms are de- 
posited onto the surface of a simple cubic crystal (i.e. 
a square lattice in d = 1) with a rate of F atoms per 
unit time and area. Atoms with no lateral neighbors are 
allowed to diffuse with diffusion constant D. In the sim- 
plest variant of the model atoms with lateral neighbors 
are assumed to be immobile so that, e.g., dimers are im- 
mobile and stable. Therefore, i* = 1. Extension of the 
rules from i* = 1 to i* = 2 is straightforward: Dimers 
are allowed to split into two adatoms, while trimers are 
stable. In either case, growth commences with a flat 
substrate, h(x, 0) = for all sites x. On deposition at x, 
h(x, t) is increased by one. We neglect barriers to inter- 
layer transport (Ehrlich-Schwoebel barriers p7|). Thus 
the only parameter of the model is the ratio D/F. This 



model was simulated in systems of size L = 65536. 

The island distance was measured here as follows. Dur- 
ing deposition, the number of islands was calculated for 
times < Ft < 1. An island was defined as a connected 
region where h(x) > a±, where a± is the vertical lattice 
constant. For early times, the number of islands increases 
due to nucleation. At roughly Ft ~ 0.5 the number of 
islands decreases again since the islands coalesce. This 
means that the island distance has a minimum which was 
used for evaluation. 
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FIG. 2. The island distance as a function of D/F for 
i* = 1 (circles) and i* = 2 (squares). The solid lines with 
slopes 7 = 0.246 and 7 = 0.288 represent fits to the last 
three decades of data. The inset shows how the local slopes 
approach the theoretically predicted exponents (dotted lines) 
for large D/F. 

In Fig. ||, the measurements of the island distance as a 
function of D/F for the MBE model with i* = 1 and i* — 
2 are shown. From the effective exponent j e g = log(10 • 
D/F) — \og(D/F) shown in the inset, it becomes clear 
that the asymptotic behaviour I ~ (D/F) 1 is reached 
only for large values of D/F. The dotted lines indicate 
the theoretical values from Eq. (||), j(i* = 1) = 1/4, and 
Eq. (||), j(i* = 2) = 2/7, respectively. Note also that the 
asymptotic exponent is approached from below. 



V. DISCUSSION 

In Fig. |, the comparison of the simulation results of 
the CGMBE model (Fig. |) and the MBE model (Fig. §) 
with various analytical formulas is shown. The small dif- 
ference between 7 measured in the CGMBE and the MBE 
model is not only due to statistical noise, but also due to 
the asymptotic behaviour of £(D/F) being reached with 
negative and positive curvature, respectively. There- 
fore, even better agreement between the results of the 
CGMBE and the MBE models for larger D/F can be 
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expected. As a consequence, the values for higher i* , 
obtained with the CGMBE model, can be regarded as 
reliable. 



length in one dimension for critical island size i* > 2 and 
have confirmed the result by computer simulations. 
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FIG. 3. The submonolayer exponent 7 as a function of i* 
in d — 1. Squares and diamonds are the results of the simula- 
tions of the CGMBE and the MBE models, respectively. The 
solid curve is the theoretical prediction for i* > 2, Eq. 
The dashed curve corresponds to formula (^) and the dotted 
line are the exponents expected, if islands up to size i* can 
diffuse, Eq. (jfij). 

The agreement of the numerical results with the the- 
oretical prediction (^|) is very good. For comparison, we 
also show Eq. (g), which in one dimension applies only 
for i* = 1. It is interesting to notice that the exponent 
for i* = 2 coincides with the one expected, if islands up 
to size i* can diffuse without decaying ||[L8) 
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VII. NOTE ADDED IN PROOF 

Following the suggestion of one of the referees we note 
that our rate equation approach applies to the general 
case of surface dimension d and island dimension d / . The 
fractal island dimension can be smaller or larger than the 
surface dimension which is integer in physically relevant 
cases. In this general case we use Ft£ d = t f with I = 
l- x l d to obtain 
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dv{i* + 1) + df : 



where v = 2 for d = 1 and v = 1 for d > 2. This 
does not apply to the exceptional case d = i* = 1 where 
d < d c = 2/i*. There equation (6) should read I = DA 2 1 
and one finds 7 = 1/(3 + df). 
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(d + 2)i* + d f 3i* + l 



(11) 



In fact, unstable islands effectively diffuse in our models. 
Consider for example i* — 2, where dimers are unstable. 
When a dimer decays, the adatoms still remain close to 
each other, and in one dimension the probability is partic- 
ularly high that they reunite again at a shifted position. 
The specific way we implemented the decay of a cluster 
in our computer simulations makes this tendency very 
clear: Whenever a dimer has formed it must break up 
in the next diffusion step, which means that each of the 
two atoms moves to one of the neighboring cells. With 
probability 1/2 they will be found in the same cell. If no 
further adatom was met there, this amounts simply to a 
displacement of the dimer by ±Ax. This analogy breaks 
down, however, for i* > 2. 

In conclusion, we have investigated the submonolayer 
epitaxial growth in a one dimensional model where is- 
lands of size up to i* are unstable while larger islands are 
stable and immobile. We have derived the exponent 7 
which determines the .D/.F-dependence of the diffusion 
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